Quintessence dynamics with two scalar fields and mixed kinetic terms 
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The dynamical properties of a model of dark energy in which two scalar fields are coupled by a 
noncanonical kinetic term are studied. We show that overall the addition of the coupling has only 
minor effects on the dynamics of the two-field system for both potentials studied, even preserving 
many of the features of the assisted quintessence scenario. The coupling of the kinetic terms en- 
larges the regions of stability of the critical points. When the potential is of an additive form, we 
find the kinetic coupling has an interesting effect on the dynamics of the fields as they approach 
the inflationary attractor, with the result that the combined equation of state of the scalar fields 
can approach —1 during the transition from a matter dominated universe to the recent period of 
acceleration. 



I. INTRODUCTION 

One goal of cosmology is to understand the origin of 
the observed accelerated expansion of the universe (see 
[l| and references therein). To date, there are several 
suggestions, including the cosmological constant, slowly 
evolving scalar fields and modifications to Einstein's the- 
ory of general relativity. Since scalar fields are predicted 
by many particle physics theories, scalar field models of 
dark energy, such as quintessence 0, d, or k-essence 
@i H 13 have been studied in considerable depth in the 
past. One requirement a satisfactory model of dark en- 
ergy must fulfill is that it leads to an equation of state 
(EOS) close to w = — 1, in order to agree with current ob- 
servational data. For single scalar fields, exponential po- 
tentials with slope A lead to a scalar field dominated uni- 
verse with late-time accelerated expansion if A < V2. In 
this case, the duration of the matter dominated epoch de- 
pends on the initial conditions for the scalar field. Thus, 
the situation is not better than that with a cosmologi- 
cal constant. If, on the other hand, A > ^3(1 + w), the 
scalar field scales with the dominant fluid (with EOS w 
d, Q). This would help the initial condition problem, 
but unfortunately, the scaling solution and the accelerat- 
ing solution are mutually exclusive. The situation with 
inverse power-law potentials is better in the sense that 
a wide range of initial conditions end up with the same 
cosmology at late times, but in order for the theory to 
be consistent with observational data, the exponent has 
to be small. 

On the other hand, if quintessence is indeed the sce- 
nario realized in nature, the quintessence sector might 
turn out to be rather nontrivial. There might, for ex- 
ample, be several scalar fields interacting p^lTll. [T^ [Tsj 
and/or the scalar fields might interact with matter in 
the universe [l^, [H, |T3|- In this paper we will 
study interacting scalar fields as a model for dark en- 



ergy and study whether this can shed some light on 
the issues in quintessence model building. Models with 
multiple scalar fields have been considered in the past: 
partly to study isocurvature (entropy) perturbations (e.g. 
ISHHIlli) and/or non-Gaussianity generated dur- 
ing inflation (see e.g. [S^, [2^ [H, HH and references 
therein). It was also found that the cumulative effect 
of many scalar fields could relax the constraints on the 
inflationary potential 26, 27, 28]. These ideas have also 
been used in models of dark energy (assisted quintessence 
[SilUjIsBl)- Like the case with inflation in the very early 
universe, there is no reason why dark energy is not driven 
by several interacting scalar fields. In a scenario with sev- 
eral scalar fields one might hope to find an explanation 
for the coincidence problem, i.e. why dark energy domi- 
nates today and not earlier during the cosmic history. In 
single scalar field models this is rather difficult to achieve 
(see e.g. fl). 

We consider a simple extension of the standard case 
with canonical normalized fields by allowing for a cross- 
term in the kinetic energy of the scalar field, which, in the 
case of a homogeneous and isotropic universe, is propor- 
tional to (/)%. This is the simplest extension possible in 
the two-field case without changing the properties of the 
potential energy. Kinetic interactions between multiple 
scalar fields have been previously studied in the context 
of models in which the dark energy equation of state can 
become less than —1 31. .32.]. 

The paper is organized as follows: in the next section 
we present the model and discuss some aspects of the 
dynamics of the system, such as effective exponents. In 
Sees, mil and HVl we perform a critical point analysis for 
two types of exponential potentials and in Sec. |V] we 
present the results of a numerical analysis. In Sec. IVII 
we consider the case of a varying coupling function. Our 
conclusions can be found in Sec. IVIII 



II. THE MODEL 
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The general action for theories which describe the dy- 
namics of two scalar fields with noncanonical kinetic 
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terms can be written as 

S = J d*Xy/^ 

where m,n = 1,2 and 



mn 7 V^m ) 



(1) 



(2) 



in units where Afpi = 1/ VSttG = 1 . In this paper, we will 
focus on a model in which a modification to the standard 
two-field system, determined by the coupling x), is 
introduced. The scalar field Lagrangian becomes 



(3) 



where, for simplicity, we have introduced the notation 
— (j) and (j)2 — X- III the following we specialize to 
a flat, homogeneous and isotropic universe, where the 
expansion is described by the scale factor R{t). 

If a is constant, one can perform a field redefinition and 
diagonalise the kinetic terms so the action takes a canon- 
ical form. This is not possible in the general case, which 
is considered in Sec. IVII We will work in the field basis 
defined above, which facilitates comparison with previ- 
ous work [l3| ■ In a similar way to the analysis presented 
in [3^ Is^l , solutions with constant a can be considered 
instantaneous critical points, valid for a particular value 
of a. 

In order to study scaling solutions in this unusual sys- 
tem we consider the scalar field system in conjunction 
with a perfect matter fiuid of energy density p, pressure 
p and EOS p = (7 — l)/5. The Friedmann equations are 



+ ^x^ + 1'i)x + v + p, 



H 



+ X +a0X + 7P), 



(4) 
(5) 



and the variation of the action with respect to the fields 

yields 



3ff 



X + 2,h{x+U]+U + V, 



0, 



(6) 
(7) 



where = dV/d(j), = dV/dx and H = R/R. The 
dot indicates a derivative with respect to time t. The ap- 
pearance of second time derivatives of both scalar fields 
in each equation is characteristic of systems of this type. 
However, in the case at hand the equations can be sepa- 
rated for ^ 4, yielding 



(l~^)+3i/<^(l-4) 
iHx[ 



X 



1 



0. 



(8) 
(9) 



To understand the behaviour when a = ±2, consider the 
reparametrization = ^(V'l +4'2), X = ^("01 ~ V'2)- Sub- 
stituting these expressions into Eqs. ^ and ([7]) gives 

(1-f ^) (V^i +3H^i' 



(1 
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FIG. 1: Assisted case parameter space showing the stability 
of the critical points for (a) a = —1 and (b) a = 1 with 7 = 1. 
The shaded area indicates the region with fJfleids = 1- When 
a is negative there are six stable critical points, defined by 
the lines 2\ + = Q, 2p + a\ = and A^g = 87. As a —> 0, 
points E and F dominate the parameter space and when a is 
positive, only these points are stable. 



It can be seen from these equations that when a = 2, ip2 
does not play a dynamical role and its value is determined 
by '(/'i. Similarly, when a = —2, ipi is not independent of 



3 







a;2 




2/2 


A 














B 


-\x2a ± \^jx\{d^ - 4) + 4 


X2 








CI 


4m 


7\/6 
2/x 







Dl 




2^6^ 





3^^(4-a2)(6-fa2^;.2) 


3{4-ai!) 






C2 


2A 


4A 


2^^(4-a2)(l-l) 





D2 


2v'6A 
3(4-a^) 


3(4-a^) 


3P^^(4 a2)(6 A2) 





E 


V 2 A 


/3"7 
V 2 M 






F 


A off 
V6A 


\/6n 


^^^(^ + |A) (1-%) 


^^A(A+|;.) (1-%) 



TABLE I: Critical points for the assisted potential 



'!/'2- Thus, when a = ±2, we are effectively dealing with 
a single field system. 

In the following we will assume that (4 - a^) > 0. 
We will perform a dynamical systems analysis, similar to 
that used in to convert the equations of motion into 
an autonomous system and investigate the properties of 
the system near the critical points. 



which the presence of additional scalar fields inhibits the 
ability of the model to produce acceleration. Defining the 
effective exponent for the soft potential by p = 2 //i^g and 
repeating the analysis above, we find 



1 _ (l-aV4) 



(15) 



A. Potentials 

We will study two forms of the scalar field poten- 
tial, the first an additive exponential potential similar 
to that used in assisted inflation [H, [S^, and assisted 
quintessence [ll,[2^[33| scenarios. 



+ e 



'MX 



(12) 



This shall be referred to henceforth as the assisted po- 
tential. The attractive feature of this potential is that in 
the case where no matter is present, one can find a solu- 
tion where the two-field system behaves like a single field 
with a potential less steep than either field. This allows 
scalar fields that are incapable of driving acceleration by 
themselves to co-operate with each other to sustain cos- 
mic acceleration. It is well known that the scale factor 
R(t) in a single field system with potential F((/)) = e~''^'^ 
evolves as i? (x V ^ where p — 2/A^. Applying a similar 
method to that used in [33] to solve for -p = 2/A^ff yields 
the effective exponent for the assisted potential in our 
system. 
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(13) 



In addition, we shall also consider a multiplicative ex- 
ponential potential, previously studied in the context of 
infiation in [HHEES, 



(14) 



This will be referred to as the soft potential, as its be- 
haviour is similar to that of soft infiation models 1381 in 



The effective exponents A^^ and can be used to sim- 
plify the untidy expressions one finds when matter is in- 
cluded. 



III. ASSISTED POTENTIAL 



After defining new variables, 



^1 -ye/f' ^2 



the autonomous system is 
+ a;iSa, 



2^2 ^a, 



y2 = -y 2^^2y2 + y2^a, 



(16) 



.2 + ^^(l-^)~'(M2/|-fA2/?) 



(17) 

(18) 
(19) 

(20) 



where ' indicates a derivative with respect to iV = In R 
and 



— —[2x1 + + 2axiX2 

+ 7(1 -xj-xl- axiX2 -yl- yl)]. (21) 
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3{4-a^) 


E 


2A + > ; 2^ + aA > ; A^g- > 87 


Stable 


37/Ae'ff 


7 


F 


2A + a/i > ; 2^ + aA > ; < A^^ < 6 


Ae'ff < 37 


1 


Apg/3 



TABLE 11: Stability of points for the assisted potential 



We define the contribution of the fields to the energy 
density of the universe by 

^f,e\ds^xl^^2+axxX2+yl + yl, (22) 

so that the constraint < fificids < 1 arises from 

The critical points corresponding to this system are 
shown in Table [H It is useful to compare these results 
to those in Table I in [l^, where the authors examine 
a system of two canonical scalar fields with the same 
potential (originally studied with nonzero curvature in 
[ll|). The points with yi, 3/2 < have been excluded as 
unphysical. It is also possible to place further constraints 
on the parameter space by excluding the solutions with 
negative values of xi and X2 i.e. where the fields roll up 
the potential. However, as there are no stable solutions 
of this kind, the unadulterated results are shown. The 
results of the analysis are summarised in Table [Til a-^d 
the stability regions of the points are plotted in Fig. [1] 

• Point A 

This is a trivial solution in which the scalar fields 
have no role to play. There are no constraints on 
the existence of the point so it is valid for all pa- 
rameter values. As rificids — and the j/i and ?/2 
values are the EOS is undefined. The eigenvalues 
are 

ei = 62 = §7 > 

63 =e4 = -|(2-7) <0, (23) 
so this point is a saddle. 

• Point B 

The point B is a kinetically driven solution with 
f^fieids = 1, given by, 

Xl = -\x2a ± \ xl{o? - A) + A 

X2 = X2, yi = t/2 = 0, (24) 



with X2 constrained only by the requirement that 
a:i,a;2 G K. The variables are real if the term un- 
der the square root is positive, giving the condition 
< (1 - a^/4)~^. Together, the -I- and - solu- 
tions form a continuous locus of points, describing 
an ellipse in the xi , X2 plane (it is a circle when 
a = 0). As ^ 4, the eccentricity increases and 
the shape rotates, anticlockwise for positive a and 
clockwise for negative a. The eigenvalues are 



ei = 3(2- 7) > 

62 =0 

63 = 3 - \J\^.X2 

64 =2,- ^\(^\x2a±\^Jl4(y^^^)^A^ .{2b) 

The ± in 64 corresponds to that in ([24]) : 61 is pos- 
itive, so the point is either an unstable node or 
saddle depending on the values of X2, a, fi and A. 

• Point CI 

The points CI and Dl have the property that only 
one of the fields' potential energy makes a contri- 
bution to the energy density. It is interesting to 
note that although the points are stable for a rela- 
tively small region of the parameter space, without 
the coupling term in the Lagrangian the points are 
always unstable. The point CI is a scaling solution 
in which the energy density of the fields decreases 
in proportion to that of the barotropic fluid. It is 
analogous to point 6 in [T5|, however, the xi value 
is nonzero and is proportional to —a. The total 
energy density of the fields goes to as — > 4. 
The condition rJficids < 1 gives the condition for 
existence, ^ > w^(4 — a^). The eigenvalues cor- 
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FIG. 2: The EOS for the assisted case with 7 = 1. (a) is a contour plot showing how Wfleids varies with jj, and A for a = 0. For 
large values of A and fi, point E is stable so lOfidds = 0. As a varies from —2 to 2, the contour lines spread out but otherwise 
do not change, (b) plots uifleids along the line fj, = X for a = — 1 (dot dashed), a = (solid) and a — 1 (dashed). In (c) uifleids is 
plotted against a for A = = 1 (solid), A = p = 1.5 (dot-dash) and A = /x = 2 (dashed). 



responding to this point are 



ei =-|(2-7)<0, 



62 
63,4 



3(2-7) 
4 



1± 



/672(4-a2) + ^2(2_9^) 
m'(2 - 7) 



(26) 



Using the condition for existence, one can show 

67'(4-a2) + M'(2-97)<A*'(2-7), 

so 63,4 < -|(2 - 7)(1 ± 1) < 0. Thus, CI is stable 
when 62 is negative i.e. /i < — aA/2. 

• Point Dl 

Dl is a scalar field-dominant solution (fificids — 1) 
analogous to 5 in so the only existence con- 
dition comes from 2/2- Therefore, the point exists 



when ^ < y|(4— a^). The eigenvalues for this 
point are 



/i(2/i + a\) 



(4 - a2) 



63,4 



2^2 



(4-a2) 



(4 - a2) 



- 3 < 



37 



(27) 



where the existence condition ensures that 63,4 are 
negative. If the point is to be stable both 61 and 
62 must be negative, so the point is a stable node 
if 



• Points C2 and D2 

These points are related to CI and D2 by the trans- 
formation xi X2, yi <-> 2/2, A ^ fi. Thus, the 
existence and stability conditions are in terms of A 
instead of fi. 

• Point E 

The point E is a scaling solution analogous to 8 in 
[l3l |. yi and 2/2 are square roots and OfioWs 7^ 1 so 
there are three conditions that have to be satisfied 
for the point to exist. To ensure that yi and ?/2 are 
real and nonzero, one has to impose 2fj, + aX > 
and 2A + a/i > 0. The third condition is given in 
terms of Aeff by X^g > 37. The eigenvalues for this 
point are 



ei.2 



63,4 



(2-7) 



1± 



2472/A2j, + 2-97 



(2-7) 



r(2-7) 



1± Wl 



87(2/x + aA)(2A + a/i) 
AA*(2-7)(4-a2) 



(28) 



(29) 



Using the existence condition, one finds + 2 



97 < (2—7). Substituting this into gives 61,2 < 
-|(2 - 7)(1 ± 1) < 0. It is clear that the term 
under the square root in (|^ must be less than 1 
if both eigenvalues are to be negative. This means 
the point is always stable when the conditions for 
existence are satisfied. 



/i < mm 



— gA 
2 



V'x(4-«') 



• Point F 

F is a scalar field-dominant solution analogous to 7 



6 



in The yi and y2 values are 



yi 



2/2 



^^/A?.A(A + |,.)(l-f 



(30) 



SO for the point to be real we require 



A,2ff(2^ + aA)(l-^ 



> 

> 0. 



Assuming ^1 ~ ~6^) < ^ ^'^"^ using (fT^]) leads to 
a contradiction, so the conditions for the point to 
exist are, 

< A^ff < 6, (2/i + a\) > 0, (2A + afi) > 0. (31) 
The eigenvalues corresponding to this point are 



62 



37 



1± 



^ 4Agff (2/i + aA)(2A + a^) 



(32) 



63^4 are negative if the term under the square root 
is less than 1, which is always true as the factor in 
the second term under the square root is positive. 
Therefore, the point is stable when A^jj < 37. 

In the case where a — 0, Xcs is smaller than A and fi. This 
means that the points E and F are similar to the scal- 
ing and scalar field-dominant solutions of a single scalar 
field [l[ with a flatter potential, which naturally makes 
an accelerating solution more favourable (cf. assisted in- 
flationary scenarios [1^). The fact that a can be both 
positive and negative complicates this slightly but the 
general principle is the same. When a > the conditions 
for existence and stability in terms of AeS are the same as 
in [Ts'l although it must be noted that the additive factor 
of a/A/i means that Acff is generally smaller for given val- 
ues of A and /i. As a ^ —2, the stability region for the 
scalar field dominated points shrinks. When a < —1 one 
can find combinations of fields that each have a small 
exponent (so considered individually they would domi- 
nate the universe) that are out of range of the rJfidds — 1 
attractors. 

The observed value of the EOS of dark energy is very 
close to ui = — 1. A problem identified by the authors 




(a) 




(b) 



FIG. 3: Soft case parameter space showing the stability of the 
critical points for (a) a = — 1 and (b) a — 1 with 7 = 1. The 
shaded area indicates the region with f2flcids ~ 1, where the 
point D in Tables IIIII and IIVI is stable. As a increases from 
—2, the shaded area grows from zero size to form a quarter 
circle when a = 0, which is then squeezed along the line X = fj, 
as a — > 2. 
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TABLE III: Critical points for the soft potential 





Existence 


Stablity 


^fields 


7ficlds 


A 


All fi,X,a 


Unstable 





Undefined 


B 


xl < (l-aV4)-i 


Unstable 


1 


2 


C 


(M?ff/6 - 1) > 


Unstable 


1 


2 


D 


Mcff < 6 


Mcff < 37 


1 


Mcff/3 


E 


Mcff > 37 


Stable 


37/Mcff 


7 



TABLE IV: Stability of points for the soft potential 



of [13] is that of reaching a sufficiently negative value 
of iffieids in a two-field system. The extra degree of free- 
dom represented by a does little to alleviate this problem 
when we consider the asymptotic values represented by 
the critical points. It can be seen in Fig. [2]that although 
a positive value of a does lead to a value of ti^geids that 
is always smaller than the a — case, the difi^erence is 
negligible in the interesting region where Wfleids ~ — 1- 
The contour plot shows that this region is only accessi- 
ble when one or more of the exponents is small, unlike 
the large values (/i, A > 1) one would expect from fun- 
damental theory. However, as we will discuss in Sec. 
IVl the kinetic coupling affects the evolution towards the 
critical point. In particular, during the transition from 
matter domination to scalar field domination, the EOS 
reaches a minimum near Wfioids = — 1 before approaching 
the asymptotic value. 

When a is negative, the region of the (^, A)-parameter 
space for which rifidds = 1 [the grey shaded region in 
Fig. [Ija)] has a similar shape to that when a is zero or 
positive but is covered by the three points Dl, D2 and F. 
^fields and the Xi, X2, yi and j/2 values are continuous on 
the boundary. In the region governed by Dl, only the x 
field has a nonzero potential and, regardless of the second 
field, its exponent fi solely determines the character of the 
critical point. D2 is similarly dependent on only one field. 
These points account for more and more of the scalar 
field-dominant (llficids = 1) region of the parameter as 
a ^ —2. The effect of this is that whenever the value of 
A is much larger than that of or vice versa the values 
of Xl, ?/2, Wfieids etc. are sensitive only to changes in the 
smaller of the two exponents. 



IV. SOFT POTENTIAL 

We define new variables for the three-dimensional soft 
system. 



Xl 



Using these, the autonomous system is 

x[ = — 3xi 

+ xiEs, 

x'2 — — 3X2 



X 



V 



(33) 



-(Axi + ixx2)y + y^s 



(34) 

(35) 
(36) 



where 



+ 7(1 



xl 



X2 



2axiX2 
— axiX2 



(37) 



We again have the constraint < ilficids < 1 from 
(|4]) , but in this case the contribution of the fields to the 
energy density of the universe is given by 



fields 



2 I 2 I I 2 

Xl + X2 + axiX2 + y 



(38) 



The critical points for this system are shown in Table 
lllli again excluding unphysical values of y. The stability 
and existence conditions are investigated below, with a 
summary of the results in Table ITVl The stability regions 
of the points are plotted in Fig. [3] 

• Points A and B 

These points are kinetically driven and are equiva- 
lent to points A and B in the assisted case. As the 
potential terms are zero in each case, the stability 
analysis is identical to that in the previous section. 

• Point C 

C is actually two points, one with the -I- sign in Xi 
and X2 and the other with the — sign. These points 
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(b) 



FIG. 4: The EOS for the soft case, (a) is a contour plot showing how lOfleWs varies with fi and A for a — and 7 = 1. For 
large values of A and n, point E is stable so lUfldds ~ 0. (b) plots Wfleids along the line /i = A for a — ~1 (dot dashed), a = 
(solid) and a — 1 (dashed). In (c) uifleids is plotted against a for X = jj. — 0.4 (solid), A = /i = 0.8 (dot-dash) and X — ji — 1.2 
(dashed) . 



are also kinetically driven, with surprisingly com- 
plicated expressions for xi and X2- The condition 
for existence is determined by the condition that 
the terms under the square roots must be positive. 
In terms of ^cS, this means that 

- 1) > 0, 

to ensure that xi and X2 are real. The only nonzero 
eigenvalue for these points is 



ei 



3(2-7)>0, 
so the points are always unstable nodes. 



(39) 



Point D 

D is a scalar field dominated solution somewhat 
similar to point F in the assisted case, as it is stable 
for values of the effective exponent less than 87. To 
ensure the y value is real and nonzero, one requires 
/i^ff < 6. The eigenvalues are 

i(Me'ff-6)<0 



ei,2 



63 



Meff 



37, (40) 



so the point is a stable node if one applies the 
stronger constraint fi^s < 37. 

• Point E 

This is a scaling solution similar to point E in the 
assisted case. The square root in y means that 
/Ltgff > 0. The density parameter of the fields has 
the value 3j/filg so the condition rificids < 1 means 
that /i^g > 87 must be satisfied if the point is to 
exist. The eigenvalues are 




(247VMgff- 



97) 



(2-7) 



(41) 



The condition for existence gives 247^//i^g < 87, 
therefore 62,3 < and this point is always stable. 

The scalar field-dominant region of the parameter 
space for the soft system is much smaller and less promis- 
ing than the assisted case. This can be understood by 
looking at the effective exponent: ficS is small only when 
/i and A are both small, so for large values of these pa- 
rameters the potential is too steep to allow acceleration. 
There are no stable rificids = 1 solutions with /i, A > 2 
and both of the exponents have to take small values if 
one is to get a solution with WfieUs ~ — 1- As in the 
assisted case, negative values of a lead to larger values 
of Wfieids and vice versa, but this effect becomes less and 
less important when fi and A are small, as can be seen in 
Figs. Hfb) and life), as near Wfidds ~ — 1 the properties 
of the system are determined by the potential terms not 
the kinetic terms. The characteristic feature of the soft 
potential is that the coupling in the potential inhibits the 
ability of the model to produce an accelerating universe. 
It can be seen from this analysis that the addition of the 
kinetic coupling exacerbates this effect. 



V. NUMERICAL ANALYSIS 

Ideally, one would like the energy density of the scalar 
fields to scale with the dominant matter component be- 
fore causing a period of acceleration. Although, as we 
have seen in the previous section, there exist scaling so- 
lutions for this system, they are of limited phenomeno- 
logical interest. This is due to the fact that the regions 
for which scaling occurs and those that give rise to ac- 
celeration are mutually exclusive and the boundaries do 
not evolve in time. Depending on the values of the con- 
stant parameters A, and a, the system either does not 
exhibit acceleration or the scalar fields quickly dominate 
the universe. There is an interesting scenario mentioned 
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Density Parameters Evolution of and y, 




FIG. 5: Numerical integration for the case a — 0, \ — 3, fi— 1 and 7=1. As in the initial conditions {(j> — l,x = 30, x ~ 
(j) — 0) are chosen so that x is initially negligible. The upper left panel shows Jlfldds (solid) and flm (dashed). The lower left 
panel shows WflcWs; the upper and lower dotted lines show the expected values at critical points E and F respectively. The 
upper right panel shows the evolution of xi (black) and yi (grey) and the lower right panel shows the evolution of X2 (black) 
and 2/2 (grey). In both of these the dotted lines are the expected values at point E and the dot-dashed lines are the expected 
values at point F. It will be seen that the dynamical variables approach the critical point values with very little oscillation. 



in [l3[ which would alleviate this problem. The authors 
describe a situation with two scalar fields in which one 
field is initially negligible compared to the other. At first 
the potential of the latter is steep so the field scales with 
the matter component. When the second field becomes 
more important, the effective exponent decreases and the 
system evolves toward a scalar field dominated critical 
point. For comparison with the kinetically coupled case, 
this is illustrated in Fig. O If one can ignore the is- 
sue of fine-tuning, this has the advantage of providing a 
mechanism for late-time scalar field dominance; however, 
as mentioned previously, the value of Wfieids is generally 
too large to compare favourably with the observational 
values of the EOS of dark energy. 

If one imposes similar initial conditions on the two- 
field system with mixed kinetic terms, the results are 
more promising. As one can see in Fig. [S] (which differs 
from Fig. [S]only in the value of a), there are two major 
differences: the energy density of the fields is negligible 



until the end of the matter dominated era and, as the 
transition occurs, the EOS of the fields drops to — 1 before 
reaching the value given in Table. |TT1 The effect becomes 
more pronounced as a — *■ 2. A similar phenomenon in 
which the EOS reaches a minimum before approaching 
the final value has been observed recently for other scalar 
field models [SO*] , however the effect is more dramatic and 
prolonged in our case. 



This behaviour can be understood by examining the 
evolution of xi, X2, yi and y2- In the a — case and 
with the X field initially negligible, the kinetic {X2) and 
potential (2/2) terms of the x field do not play a role until 
eventually j/2 increases. In the a = 1.99 case, exp(— /iXini) 
is small so ^ V^. In this case, the equations of mo- 
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Density Parameters Evolution of and y, 




FIG. 6: Numerical integration for the case a = 1.99, A = 3, pi = 1 and 7 = 1. (The quantities plotted and initial conditions are 
the same as those in Fig. [S]) One can see that while the fields are scaling, f2flcids remains small and during the transition from 
matter to scalar field dominance, lOfleWa reaches —1. 



tions for both fields [see ([8]) and ([9])] reduce to 



a"* 



(l-a2/4) 
(l-a2/4) 



Integrating these equations, with the initial conditions 
= X = 0, we find x — ~4' ^^id so X2 ~ —xi, as can be 
seen in Fig. [6l The factor of (1 — 0^/4) in the denom- 
inators of the equations above increases the values of cp 
and X compared to the a = case so when the system 
approaches the scalar field-dominant solution they over- 
shoot the critical point and oscillate around it. As the 
term involving the coupling {axiX2) is significantly large 
and negative, Wficids is decreased during this period. 

Assuming X2 — —xi, yi, j/2 <C 1 and a « 2 means that 
the function Ea (defined in pT|) ) simplifies to Sa 
This means that (fT9)l becomes 



§7- 



Vi 



2 A 

37 



Xi 



(42) 



As in the a = case, the system first approaches the 
scaling point E, which means xi ^f^l/^- Thus, the 
evolution of yi is quickly halted and, until the x-ficld 
becomes significant, the fields only occupy a fraction of 
the energy density of the universe. Effects of this type are 
not seen in the case of the soft potential as the condition 
Vx <C Vif, corresponds to /i 3> A, which is not realised in 
the region where the rJficids = 1 attractor is stable. 

This feedback effect occurs whenever a is close to 2 
and ^ (or <^ V^) so is affected both by the 
initial conditions and the relative sizes of A and /i. It is 
noteworthy that one does not actually require the initial 
conditions of the two fields to be different to see the scal- 
ing effect and the decrease in Wfidds (see Fig. [T]). If the 
initial contribution of the scalar fields to energy density 
is small compared to that of matter, a small difference in 
the values of jj, and A will yield the condition xi « —X2 
that pushes the EOS to —1. When the system is initially 
close to the critical point there is no scaling regime, but 
as long as ^ and a w 2, one still observes the effect 
on Wficids- However, although the scalar field dominated 
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(a) (b) 

FIG. 7: Two trajectories for different initial conditions {(ph-ii ~ Xini), for X — 3, fi = 1 and 7 = 1. The upper plots show the 
evolution of the energy densities of the fields (solid) and matter (dashed) and the lower plots show the scale factor, luflcids. 



solution is an attractor, the question of when the scalar 
fields become cosmologically important is dependent the 
initial conditions. 

We can gain further insight into this behaviour (in the 
case of constant a) by performing a field redefinition to 
diagonalise the kinetic terms, 



^1 



f2 



1 



f) 



1+f) 



x), 
x), 



(43) 
(44) 



so that the Lagrangian for the scalar field system with 
the assisted potential becomes, 



C 



2 V^MV 



(45) 
2, 



with A = l/y/2 - a and B = 1/^2 + a. As a 
B ^ 1/2 and ipi ^ [note that the combination Aipi — 
i(</) — x) remains finite]. Thus the Lagrangian can be 
written 



C 



■i {d^V2d''ip2) - - T/ae-'^^^/^^ (46) 



with Vi = e-^^f^ and V2 = e'^^'^i constant. This is 
similar to the case in [lo| . where a double exponential 
potential is used to give a scaling regime followed by ac- 
celeration driven by the scalar field. 



VI. FIELD DEPENDENT COUPLING 



In the previous analysis, we assumed that the kinetic 
coupling, a, is constant in time and independent of the 
scalar fields. In this section we relax this assumption 
and extend the analysis to include a = a{<j)^x) for the 
assisted potential. The Friedmann equations Q and ^ 
are unaffected but the equations of motion ([5]) and ^ 
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FIG. 8: A slowly varying coupling function, modelled by 
a(0,x) = logio (|(<^ + X) + 5) - ^, with cj!>i„i = Xmi, A = 2, 
fj, = 2.1 and 7 = 1. The upper panel shows fifleids (solid) and 
Qm (dashed) and the middle panel shows the EOS (solid) 
and the value of Wfidds at the critical point, (wfioids)F (dot- 
ted). The lower panel shows the evolution of a. The system 
initially heads toward the field-dominant point F. As a in- 
creases, lOfleids tracks the change in (uifleids)F. 

have additional terms on the right-hand side, 
4> + 3H<j> + (l - s^y\v^ - ^V^) 

= (l-^)"'(W0'-iaxX'), (47) 
X + 3HX + (l-x) '(^x-fl^0) 

= (l-x) '(i««xX'- W), (48) 

where and are the partial derivatives of a with re- 
spect to the fields <j) and x respectively. Our previous 
analysis is valid for two cases of particular interest: a 
slowly varying coupling and a rapid change between ap- 
proximately constant values of a. In the former case, the 
effect of the extra terms on the dynamics of the fields is 
very slight and we can consider the asymptotic solutions 
discussed in Sec. IIIII as instantaneous critical points (cf. 
dEllll)) vahd for particular values of a((/), x). The sys- 
tem will evolve toward one of the stable critical points 
as before; as a varies only slowly, the solution will track 
the variation in a. This can be seen in Fig. [51 where the 
solution tracks the point F as a slowly increases to larger 




N 



FIG. 9: A rapid shift in a, modelled by a{(j>,x) ~ 1-9 x 
ta.nh{8[{(f) + x) - 25]), with (/)i„i = Xmi, A = 2, = 2.1 and 
7 = 1. The quantities plotted are the same as those in Fig. [8] 
except that when point F is unstable, the dotted line in the 
second panel shows (u'ficids)£; ~ 0. Initially the point E is sta- 
ble and the system exhibits scaling. After a changes rapidly 
the existence conditions for E are no longer satisfied and the 
system evolves toward the field dominant point F. 



values. 

When the rate of change of a is faster than the typical 
response time of the system, in order to understand the 
behaviour one would typically have to study the effect 
of the terms in and in (|47|l and ((48|) for a par- 
ticular model. An interesting exception to this is when 
the timescale during which a changes is short. An exam- 
ple of this is illustrated in Fig. [9l where we have chosen 
a{(j),x) oc tanh(C(0 + x) + constant). This allows us 
to generate a rapid change in the coupling a. At first, 
a is approximately —1.9 so the system exhibits scaling. 
When the system has reached the scaling point (E in this 
case), Wficids = and there can be no acceleration. The 
scaling regime is abruptly ended when a — > -1-1.9. Point 
E no longer exists and F is stable instead, so the scalar 
fields will drive the acceleration of the universe. 

To conclude, the results of the phase space analysis 
in the case of constant a are useful when a is a slowly 
evolving function as well as the case of a spontaneous 
change in the value of a (followed by a period where the 
coupling is approximately constant). 
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VII. DISCUSSION 

Although the use of scalar fields to model the acceler- 
ation of the universe could well be described as standard 
practice in cosmology, the fundamental physics that un- 
derpins these phenomenological models is by no means 
clear. It is wise, therefore, to consider the possibility that 
the observed acceleration of the universe may have more 
exotic origins. The recent focus on scalar fields with non- 
canonical kinetic terms derived from string-inspired mod- 
els has shown that models that deviate from the standard 
lore are capable of producing fruitful results that display 
qualitatively different behaviour. In this work we have 
considered a model of dark energy in which the kinetic 
terms of two scalar fields are coupled by a term in the 
Lagrangian of the form a{da(j)d°'x)- This allowed us to 
study the simplest extension of the standard case with- 
out modifying the potential energy. After analysing the 
phase plane structure of the model for two different po- 
tentials, we found that the basic features of the uncou- 
pled case were preserved and the most important effect 
of the extra degree of freedom was to change the range of 
parameters that lead to the different types of solutions. 
In particular, there exist stable points corresponding to 
scalar field-dominant or scaling solutions for every value 
oi X, fi and a considered. 

In the case of the assisted potential we found that for 
negative a, critical points CI, Dl, C2 and D2 that existed 
but were unstable in the case without the kinetic coupling 
became stable. This means that when a goes from pos- 
itive to negative values, there are three different regions 
in the A, /i-parameter space that give rise to accelerated 
solutions, each with a different value of the EOS, Wfioids- 
This is less interesting from a phenomenological point of 
view, however, as except when the values of the potential 
exponents are very small, negative a values have Wficids 
larger than the observed dark energy EOS. In the case 
of positive a, one finds the EOS to be closer to —1 than 
that in the model with a = 0, though in the region of 
interest (wfioids ~ —1) the effect of the kinetic coupling 
is extremely slight. The soft potential evinced a similar 
dependence on a, though it is of less interest as a cosmo- 
logical model as accelerated behaviour only occurs with 
small values of the exponents in all cases. 

Although models with exponential potentials exhibit 
both scaling and scalar field-dominant solutions, it is 



difficult to solve the coincidence problem without fine- 
tuning the initial conditions or allowing the parameters 
to vary in time. In Sec. |V] we discuss how the kinetic 
coupling affects the dynamics of the fields as they evolve 
towards the critical points. When a is close to 2 one finds 
that the value of Wficids becomes very close to —1 for a 
brief period before the scalar fields dominate. This is a 
partial resolution of the perennial problem of the EOS 
being too large. 

This analysis is revealing as the addition of a cross- 
kinetic term has only a relatively minor effect on the 
dynamics of the two-field system for both cases stud- 
ied, even preserving many of the features of the assisted 
quintessence scenario. In the case of the assisted poten- 
tial, the kinetic coupling has an interesting effect on the 
dynamics of the fields as they approach the stable so- 
lution, with the result that the EOS of the scalar fields 
can approach —1 during the transition from a matter 
dominated universe to the recent period of acceleration, 
something that does not occur in assisted quintessence 
or the single field case. A natural extension to this work 
is to consider the case when the scalar fields are coupled 
directly to the matter fluid by a similar mechanism to 
that described in [l^. In this case, one would expect the 
value and stability ranges of the critical points to change 
dramatically due to the extra degree(s) of freedom, but it 
would be interesting to see whether the additive poten- 
tial still allows large values of the potential exponents. 
The work could also be extended by treating A and /i as 
dynamical variables in a manner similar to |34| , in which 
case it may be possible to obtain solutions with a large 
range of initial conditions that track the background so- 
lution and later dominate. 

As we have discussed, our results are also applicable in 
the case of a slowly varying coupling function a{(j),x) 
well as the case in which the coupling undergoes a rapid 
change but is constant (or evolving very slowly) after the 
transition. 
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